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Abstract. Searches for gravitational waves with km-scale laser interferometers 
often involve the long- wavelength approximation to describe the detector 
response. The prevailing assumption is that the corrections to the detector 
response due to its finite size are small and the errors due to the long- wavelength 
approximation are negligible. Recently, however, Baskaran and Grishchuk (2004 
Class. Quantum Grav. 21 4041) found that in a simple Michelson interferometer 
such errors can be as large as 10 percent. For more accurate analysis, these 
authors proposed to use a linear-frequency correction to the long-wavelength 
approximation. In this paper we revisit these calculations. We show that 
the linear-frequency correction is inadequate for certain locations in the sky 
and therefore accurate analysis requires taking into account the exact formula, 
commonly derived from the photon round-trip propagation time. Also, we 
extend the calculations to include the effect of Fabry-Perot resonators in the 
interferometer arms. Here we show that a simple approximation which combines 
the long- wave length Michelson response with the single-pole approximation to the 
Fabry-Perot transfer function produces rather accurate results. In particular, the 
difference between the exact and the approximate formulae is at most 2-3 percent 
for those locations in the sky where the detector response is greater than half of 
its maximum value. We analyse the impact of such errors on detection sensitivity 
and parameter estimation in searches for periodic gravitational waves emitted 
by a known pulsar, and in searches for an isotropic stochastic gravitational- wave 
background. At frequencies up to 1 kHz, the effect of such errors is at most 1-2 
percent. For higher frequencies, or if more accuracy is required, one should use 
the exact formula for the detector response. 



PACS numbers: 04.80. Nn, 95.55.Ym 
1. Introduction 

Searches for gravitational waves are currently conducted with km-scale laser 
interferometers such as LIGO [1 and VIRGO [2]. These detectors utilize a Michelson 
configuration which is further enhanced by the addition of Fabry-Perot cavities in 
the interferometer arms. Development of efficient data analysis algorithms requires 
accurate knowledge of the response of these detectors to gravitational waves. There 
are two somewhat different points of view on how to calculate the detector response. 
In one approach, it is assumed that the size of the detector is much less than the 
wavelength of the incoming gravitational wave, and therefore can be neglected in the 
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calculations. This approach is often called the long-wavelength approximation [314]. 
The advantage of this approximation is that it allows one to interpret the effect of 
gravitational waves entirely in terms of the motion of test masses, which is appealing 
to our physical intuition. In another approach, one takes into account the finite size 
of the detector by considering variations in the gravitational wave within the duration 
of one photon round trip between the test masses [5 -[10 . The detector response 
obtained in this way is more accurate but no longer allows the simple interpretation 
in terms of test mass motion [TTj. Such calculations are commonly used to derive 
the response of space-borne gravitational- wave antennae [T?-p!6^. For ground-based 
detectors, one usually adopts the long- wavelength approximation, assuming that it is 
sufficiently accurate. 

It was pointed out by Baskaran and Grishchuk p/7| that even for ground-based 
detectors, the long- wavelength approximation can lead to noticeable errors in the 
estimation of parameters of a gravitational wave. In particular, they found that 
the error in searches for periodic gravitational waves can be as large as 10%, thus 
raising a concern about the validity of the long-wavelength approximation in recent 
searches for gravitational waves. Some key points of their analysis however required 
clarification. The authors assumed that for ground-based detectors it suffices to use 
the first-order correction to the long- wavelength approximation and thus introduced a 
linear- frequency detector response, whereas the exact, non-linear formula was readily 
available. Also, the calculations did not take into account the presence of Fabry- Perot 
cavities in the interferometer arms, which play a crucial role in the formation of the 
signal. It is therefore worthwhile to reconsider this analysis. 

In this paper we re-evaluate the errors due to the long- wavelength approximation 
and assess their impact on current searches for gravitational waves with km-scale laser 
interferometers. We show that the linear- frequency approximation is inadequate for 
some locations in the sky, and therefore one must use the exact formula for the detector 
response to estimate systematic errors from the long- wavelength approximation. To 
make the analysis applicable for LIGO and VIRGO detectors, we include the transfer 
function of Fabry- Perot cavities in the interferometer arms. Using the exact expression 
for the detector response, we estimate the errors resulting from the long-wavelength 
approximation in searches for periodic gravitational waves and in searches for an 
isotropic stochastic gravitational-wave background. 



2. Michelson interferometer response (long-wavelength approximation) 

In the transverse-traceless gauge ^ , a plane gravitational wave coming from direction 
h on the sky is given by 

hij{t,x) = /i+(t,f)e+.(n) ^hy,{t,x)e^j{h) , (1) 

where /i+^x {t^ x) = hj^^^{t + x • n/c), and the polarisation tensors are 

= ^i^j - mrrij, (2) 

The unit vectors £ and rh are chosen so that m, h form a right-handed orthonormal 
basis. The rotational degree of freedom associated with the choice of £ and rh in the 
plane perpendicular to n is often called the polarisation angle p. In what follows, we 
will suppress /^-dependence for simplicity. 
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Consider a Michelson interferometer with arms aligned along the unit vectors a 
and b. In the long-wavelength approximation, a signal produced by a gravitational 
wave in the detector [H [18] is given by 

F(t) = i(aV-6V)%(t,0), (4) 

where we assumed that the detector is located at x = and its size is negligible. 
Equivalently, the signal can be written as 

V{t) = F^{n)h^{t) + Fx(n)/ix (t) , (5) 

where 

FA{n) = ^ {aiaj - hihj)e'l{n) (6) 

are the interferometer responses to the two independent polarisations (A = +, x) of 
the gravitational wave. In the frequency domain, ([5| becomes 

VU)=F+{h)h+{f)+Fy,{n)hy.U). (7) 

In what follows, tilde always denotes Fourier transform with respect to t. 

Three-dimensional representations of the absolute value of Fa as a function 
of n are often called antenna patterns. Antenna patterns have traditionally been 
shown for a particular choice of polarisation basis: £ = and rh = (j)^ where 6 
and (/) are the unit vectors corresponding to the spherical coordinates (j) G [0,360°] 
and e [0,180°]. Figure [l] shows the antenna patterns in the coordinate system 
with the X and y axes aligned with the interferometer arms. In these coordinates, 
h = (sin cos (/), sin sin (/), cos 0) . 

3. Michelson interferometer response (exact formula) 

The detector response which takes into account the finite size of the interferometer 
will be called here exact in contrast to the approximate response (|6|. Here we give a 
brief derivation of the exact detector response following recent calculations in [19l [20] . 
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3.1. Photon propagation time 

The interval for photons propagating in spacetime with a gravitational wave ([T]) is 

= -c^dt'^ + [Sij + hij {t, x)] dx'dx^ = . (8) 

Consider a photon launched in the direction d to be bounced back by a mirror some 
distance L away. On the way forward, the unperturbed photon trajectory is = a*^, 
where ^ G [0,L]. Substituting this trajectory in ([8| and solving for t, we obtain 

c{t - to) = (1 + hij a'a^)^'^ d^! . (9) 

^0 

Let T be the nominal (unperturbed) photon transit time: T = L/c. In the 
presence of a gravitational wave, the transit time will slightly deviate from its nominal 
value giving rise to a small perturbation: 

6T{t) = ^^a'a^ hij + ^ + d^, (10) 

where to is the starting time for the photon propagation which can be approximated 
hy to = t — T. Similarly, on the way back, 

ST'it) = ^^a'ai % [to + ^ + d^, (11) 

where to can also be approximated by to = t — T. Then the perturbation of the 
round-trip time is given by 

ST,,tXt) = 5T(t - T) + 5T'(t). (12) 

In the Fourier domain, it can be written as 

= a^aj D{a, /) 4{n) IaU), (13) 

where we introduced the transfer function 

D{a, /) = [e^^/^+sinc (7r/T_) + e-'^/^-sinc (7r/T+)] , (14) 

with short-hand notation: T± = T(l ± a • n). Further calculations require switching 
from photons to continuous electro-magnetic waves. 



3.2. Phase lag of a continuous electro -magnetic wave 

For an electro-magnetic wave propagating in the a-direction, the electric field is given 
by £{t^x) = Aex.-p[i{u;t — kx ■ d)]^ where A is the amplitude, co is the frequency, and k 
is the wavenumber {k = co/c). It is convenient to suppress the fast-oscillating factor 
^tujt introducing the slowly- varying amplitude [21 : E = fe~*^^. Consider a simple 
Michelson interferometer with equal arm lengths, L, as shown in figure [2] (left). The 
phase delay of the electro-magnetic field returning to the beam splitter after a round 
trip in the arm is a;(2T + STj^^.) = + ?/^, where 

^{t)=uST,,tXt) (15) 

is the phase delay due to the gravitational wave. Two such phases corresponding 
to arms d and b will be denoted here ipa and ip^^. Let the amplitude of the field 
immediately after the beam splitter be Ao. Then the amplitude of the field incident 
on the beam splitter from arm a is = ^0 e^p{—2ikL — iipa), and similarly for E^. 
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Figure 2. Simple Michelson interferometer (left) and Michelson interferometer 
with Fabry-Perot arm cavities (right). Enlarged is a schematic picture of the 
photon trajectory in the d arm. (Forward and return paths are separated for 
clarity.) Here we neglected the recycling mirror which increases the power incident 
on the beam splitter but otherwise does not affect the detector response. 



v{t) = —[Mt)-Mt)]- (17) 



Assuming that the interferometer operates at the dark fringe (destructive interference), 
we find that the field at the output (signal) port is proportional to 

Ea{t) - E-,{t) « -iAoe-^^'^^lMt) - Mt)] ■ (16) 
With appropriate normalizatior||] the signal is given by 

1 

In the Fourier domain, it can be written as 

V{f) = G+{h, f)h+{f) + Gx (n, f)h^ if) , (18) 

where Ga{^^ f) are the exact detector responses to the two independent polarisations 
of the gravitational wave: 

GA{n, f) = \ [ata, D{a, /) - hb, D{b, /)] e^(n) . (19) 

Note that the long- wavelength formula (|6| is a special case of the exact response: 

FA{n) = GA{n,0). (20) 

As the frequency of the gravitational wave increases, the difference between Fa and 
Ga becomes more and more pronounced [9l [22] . The most drastic change occurs at 
the inverse of the photon round- trip time: / = 1/(2T). In Fabry-Perot cavities this 
quantity is called the free spectral range or FSR (see section [XT]) . 

Figure [3] shows the magnitude of the detector response functions Ga(^,/) at 
the free spectral range of the 4-km LIGO interferometers (/ = 37.5 kHz). Note the 
presence of additional lobes in the antenna patterns, and a factor of 5-8 reduction in 
magnitude compared to the response at / = 0, shown in figure [T] 

X The normalization is such that max^(GA) = 1 at / = for both A = +, x. 
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Figure 3. Antenna patterns at the FSR frequency (37.5 kHz): |G+(ri, /)| (left) 
and |Gx(ri,/)| (right). 



3. 3. Approximate formulae for the Michelson response 

The long-wavelength approximation (|6| is obtained by entirely neglecting the 
frequency dependence of the exact response 

GA{n,f)-FA{h). (21) 

A moderate frequency dependence is obtained by adding the first-order correction: 

GA{n,f)^FA{h) + f^{n,0). (22) 

This may not necessarily be a better approximation than ([21]) since the higher order 



terms neglected in (22) can be greater than the first-order term (proportional to first 
derivative). This happens because the first-order term vanishes at some locations in 
the sky whereas the sum of all higher order terms remains non-zero. 

The most recent estimation of errors introduced by the long-wavelength 
approximation is due to Baskaran and Grishchuk ^Jj. These authors argued in 
favour of the linear approximation of the detector response ( [22| , and estimated 
the corrections to the long-wavelength formula using the first-order (/-proportional) 



term. (In [17], the two terms in the right-hand side of (22) are called electric and 
magnetic components.) As we have shown, such a linear approximation is problematic. 
Fortunately, one need not be concerned with the accuracy of the linear approximation 
since the exact formula is readily available. 

4. Transfer function of Fabry-Perot arm cavities 

The response of LIGO and VIRGO detectors is enhanced by incorporation of optical 
resonators (Fabry-Perot cavities) in interferometer arms. Here we briefly derive the 
transfer function of a Fabry-Perot cavity and include it in the exact detector response. 

4.1. Phase amplification due to multi-beam interference 

Consider a Fabry-Perot cavity in one of the arms of the interferometer, as shown in 
figure [2] (right). Let the amplitude of the light incident on the cavity be Aq. Then the 
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fields circulating in the cavity satisfy 

Ei{t)=tiAo-riE2{t), (23) 

E2{t) = -r2Ei{t- 2T)e-^'^^-''^^^\ (24) 

where ti is the transmissivity of the front mirror, and ri^2 are the reflectivities of the 
front and back mirror, respectively. Here i/j is the phase lag due to the gravitational 



wave, (15). The condition for resonance implies that L is equal to an integer number 
of half- wavelengths of light and therefore e~^*^^ = 1. 

The fleld returning to the beam splitter from the cavity, = riAo -\- tiE2, 
consists of the promptly reflected fleld and the leakage fleld. The information about 
the gravitational wave is contained in the leakage fleld which is proportional to internal 
fleld E2. Let the amplitude and the phase of this fleld be A and ^, i.e. E2 = Ae~*^. 



Solving equations (|23[) and (24) to flrst order in ^, we flnd that the amplitude is given 



hy A = — tir2Ao/(l — rir2) and that the phase satisfles the equation 

^{t) - rir2*(i - 2T) = ^ , (25) 

or equivalently, 

oo 

= ^(rirs) V(^ - 2/cT) . (26) 

k=0 

Taking the Fourier transform of either (l25l) or (|26|), we obtain 



*(/) = goCimf) , C{f) = ^ _ , (27) 

where = (1 — ^1^2)"^ is the cavity ampliflcation factor and C{f) is the normalized 
transfer function. Note that C{f) is a periodic function of frequency with the period 
known as the free spectral range, FSR = 1/(2T). The 4-km LIGO interferometers 
have the cavity gain of 70.6 and the FSR of 37.5 kHz. 

It is convenient to represent C{f) in the following equivalent form: 

r(f\- J^^fT sinh(27r/oT) 

sinh[2^/oT(l + z///o)]' ^''^ 

where /o is the lowest order pole, /o = | ln(rir2)|/(47rT). At low frequencies 
(/ <C FSR), one can approximate the exact response ([28]) with a zero-pole fllter, 

where /i = FSR/tt is the frequency of the zero. In the 4-km LIGO interferometers, 
/o = 85.1 Hz and /i = 11.9 kHz. 

4.2. The response of a Michels on- Fabry -Perot interferometer 



Calculating the fleld at the output (signal) port, as we did in section 3.2, we flnd that 



the signal in a Michelson interferometer with Fabry-Perot arm cavities is 

V{f) = H+{n, f)h+{f ) + ffx (n, /)^x (/), (30) 



where Ha can be found by combining the Michelson response (19) with the transfer 



function of a Fabry-Perot cavity (27): 



HA{nJ)=CU)GA{nJ). (31) 
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Figure 4. Top row: the normalized detector responses to + and x polarisations, 
bottom row: the magnitude of the error. All these quantities are calculated at 
1.2 kHz for comparison with |17) . 



Here we omitted the cavity gain (^fo) to have the detector response normahzed to 1 at 

/ = o. 

An ad-hoc low-frequency approximation for this formula is obtained by replacing 
the exact Michelson response, Ga(^, /), with its long- wavelength counterpart, ^^(n), 
and by replacing the exact Fabry- Perot response, with the single-pole transfer 

function. 

The result is the long- wavelength approximation for a Michelson-Fabry-Perot 
interferometer: 

i^A"(n,/) = Cpoie(/)FA(n), (33) 
which is frequently used in theoretical calculations and data- analysis algorithms. 



The difference between the exact (31) and approximate (33) detector response. 



SHa^ is a source of systematic errors. The magnitude of SHa at a given frequency 
is a function of sky location, as shown in figure |4] Note that the relative error, 
ca = \5Ha\/\Ha\ diverges at those places in the sky where Ha vanishes and SHa 
remains non-zero. Since the signal also vanishes at these locations, we can safely 
exclude all such places from the error estimation. Therefore, we consider only those 
locations in the sky where \Ha\ is greater than a certain value (threshold). For a 
conservative threshold of 25% of the maximum value, we find that ca is at most 6-7%. 
With a more realistic threshold of 50%, ca is at most 2-3%. 



4.3. Variations of the interferometer arm lengths 



In the following analysis, we will also need to consider the response of the 
interferometer to changes of its arm lengths. Such a detector response is commonly 
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used in calibration measurements. If xi and X2 are the displacements of the front 
and back mirrors of one of the arm cavities, the change in the distance between the 
mirrors, as seen by the light propagating in the cavity, is 

dL{t)=X2{t-T)-xi{t), (34) 

where T is the photon transit time {T = L/c). In this case, the signaj§]is 

V{f) = f,C{f)5L{f) 

= MC(/)[e-'2-/^X2(/)-xi(/)], (35) 



where C(/) is given by (|27|. (For derivation of Q and ([35]), see [21,.) At low 
frequencies, this response can be approximated by 

V{f) « M [CpoUf] Hf) - CM • (36) 



The relative error between the exact ([35| and approximate formula (36) is less than 
1% for the front mirror and 0.5% for back mirror for frequencies up to 2 kHz. 



5. Effect on searches for periodic gravitational waves 



A gravitational- wave signal from a pulsar is quasi-periodic and therefore greatly 
benefits from synchronous detection (heterodyne method) |23j. This method removes 
the dominant oscillatory part of the signal at frequency /, and corrects for any phase 
modulation (Doppler shift) due to the rotation of the Earth and its orbital motion 
relative to the pulsar, as well as any possible pulsar frequency evolution (e.g., spin- 
down). If we ignore variations in the detector response due to these small frequency 
shifts, the output of the heterodyne method is given by the following (complex) signal 



yit) 



- cos^ l) — iHx (^, /; t) cos l 



hoe" 



(37) 



where ho is the amplitude, i the inclination angle, and <po the initial phase of the 
heterodyne transformation. There is also an implicit dependence on the polarisation 
angle /3 which enters the signal through the detector response functions. The remaining 
time dependence in (37) comes from the sidereal rotation of unit vectors a and b which 



are fixed to the Earth. 

The use of the long-wavelength approximation affects both detection sensitivity 
and parameter estimation. For simplicity, assume that all parameters of the pulsar 
signal are known except for the amplitude ho. Then the relevant quantities are 

^^'"^ and S-^^^'^ 



1 



1. 



(38) 



where y{t) is given by ([37[) and z{t) is calculated with the same expression except that 



Ha are replaced with ^^p^^^ given in (|33|). The inner product is defined as 



(^n, Vv) = J^^ [yu{t)yl{t)] dt , 



(39) 



where the integration is taken over one sidereal day. It can be shown that e is the 
fractional change in signal-to-noise ratio, and 5 is the fractional bias in the estimate 
of the amplitude of the gravitational wave, both caused by the use of the inaccurate 
detector response [24] . 



§ Consistency with a gravitational- wave signal requires that = 2/L. 
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Figure 5. e and 5 as a function of frequency for dec = 45° (left), and as a 
function of declination for f = 2 kHz (right). (Other parameters are /3 = 0, l = 0, 
and (f)o = 0.) 

Figure |5] shows e and (5 as a function of the frequency of the gravitational wave 
for fixed source dechnation, and also as a function of the declination angle for fixed 
frequency, for the 4-km LIGO Hanford interferometer. For source location, we only 
need to specify the declination angle as the dependence on the right ascention is 
removed by the 24-hour sidereal-time integration. For simplicity, we assumed that 
P = L = (j)o = 0, which corresponds to a circularly-polarised gravitational wave. 
One can see from the figure that the change in the snr, e, is much less than 1% for 
all directions on the sky and all frequencies up to 4 kHz. Also, the bias, (5, in the 
estimation of ho is less than 1.5% for all directions on the sky at / = 2 kHz. 

6. Effect on searches for stochastic gravitational waves 

Consider an isotropic, unpolarised stochastic gravitational- wave background and 
assume that it is described by a stationary Gaussian random process. Then the 
expectation value of the cross-correlation [8l ES] of the outputs of two detectors is 
proportional to the overlap reduction function: 

m = ^ [ d'^^HMn. fmAn. /)e^2./n.(.^-.-.)/c ^ (^q) 

OTT Js'2 

where xi and X2 are the locations of the two detectors on Earth and HiA{n^ f) and 
^2a(^,/) are their response functions. (The summation over polarisation indices A 
is understood.) Figure [g] shows a typical r(/) calculated with both the exact ( [31] ) and 
approximate ([33| formulae. The two detectors are the 4-km LIGO Hanford (HI) and 
Livingston (LI) interferometers. 

The error in the detector response from the long- wavelength approximation affects 
detection sensitivity and parameter estimation via the overlap reduction function. The 
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Figure 6. Overlap reduction function for LIGO HI and LI interferometers, 
calculated using the long- wavelength and exact detector response functions. Left: 
Re(r(/)) for frequency range from 1 Hz to 1 kHz. Right: a small section of the 
plot magnified to show the slight difference between the two curves. 



two relevant quantities are 



^^'^^■"'•^ and 5=-J£l%L_i, (41) 



where r(/) is given by (40) and r£.^.(/) is calculated with the same expression except 
that Ha are replaced with j^^pp^^^ given in (33). The inner product is defined as 

,p P . _ ReK(/)r„(/)] 

where (5i,2(/) are the power spectra of the outputs of the two detectors. It can be 
shown that e is the fractional change in signal-to- noise ratio, and S is the fractional 
bias in the cross-correlation statistic, both caused by the use of the inaccurate detector 
response [24 . In this calculation, we also took into account the systematic error from 
the single-pole approximation to the cavity response [see (36)] which occurs when the 
detector output is calibrated. 

Table [l] shows e and S corresponding to several detector cross correlations. The 
upper part of the table gives the values for the LIGO- ALLEGRO search for a stochastic 
background [28 , which involved correlations of the Livingston interferometer with the 
ALLEGRO bar detector in a narrow band of frequencies near its peak sensitivity at 
915 Hz. The search was performed with three different orientations of the bar detector: 
parallel to the X-arm of the interferometer (AX), parallel to the Y-arm (AY), and 
parallel to the bisector of the two arms (AN), also known as the null orientation. 
We find that the upper limits in [28] are not affected to the stated precision by the 
corrections in F. 

The rest of the table gives the values for e and 5 corresponding to various choices 
of cross-correlation between the LIGO and VIRGO (VI) interferometers. The first 
frequency band, 50-150 Hz, corresponds to the best sensitivity of the LIGO detectors. 
The second frequency band, 900-1000 Hz, is motivated by the VIRGO detector. At 
the present, this is where it contributes the most to the correlation-based searches. 
Note that only the Hl-Ll low frequency analysis has been done so far, see e.g., [29 . 
The other interferometer cross-correlations will be performed in the future. One 
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Table 1. The fractional snr change, e, and the fractional bias, S, for an isotropic 
stochastic background. Here we assumed that the stochastic background has 
constant energy density in the frequency band of interest [/min, /max], and used 
the nominal design sensitivities of the instruments. (For details, see 24^.) 
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can see from the table that ah these errors are less than 1%, except for the LIGO- 
VIRGO correlation searches around 1 kHz, which would have a 1-2% fractional bias 
in parameter estimation. 

7. Summary 

We re-evaluated high-frequency corrections to the detector response and analysed their 
effect on current searches for gravitational waves with km-scale laser interferometers. 
Using the exact formula for the detector response, we estimated systematic errors 
introduced by the long-wavelength approximation in detection sensitivity and 
parameter estimation of gravitational waves. Typical examples were taken from 
searches for periodic gravitational waves and from searches for an isotropic stochastic 
background. So far, in all cases considered, the errors are at most 1-2% and somewhat 
smaller than the previously reported 10% error at 1.2 kHz [T7] . 

We have thus shown that the long- wavelength approximation for IVIichelson- 



Fabry- Perot interferometer (33) was sufficiently accurate for searches performed to 
date, which were limited to frequencies below 2 kHz. However, extending the analysis 
to higher frequencies will likely require using the exact formula for the detector 



response (31). For example, the exact formula is essential in searches for burst [30] 
and stochastic gravitational waves [31 at the free-spectral-range frequency (37.5 kHz) 
of LIGO interferometers. Future searches for stochastic gravitational waves can go 
to even higher frequencies [32 . In conclusion, we recommend the use of the exact 
formula whenever the accuracy of the detector response is important. 
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